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Abstract. This work establishes nonlinear orbital asymptotic stability of scalar radiative shock 
profiles, namely, traveling wave solutions to the simplified model system of radiating gas [8], consisting 
of a scalar conservation law coupled with an elliptic equation for the radiation flux. The method is 
based on the derivation of pointwise Green function bounds and description of the linearized solution 
operator. A new feature in the present analysis is the construction of the resolvent kernel for the case 
of an eigenvalue system of equations of degenerate type. Nonlinear stability then follows in standard 
fashion by linear estimates derived from these pointwise bounds, combined with nonlinear-damping 
type energy estimates. 
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1. Introduction. The one-dimensional motion of a radiating gas (due to high- 
temperature effects) can be modeled by the compressible Euler equations coupled with 
an elliptic equation for the radiative flux term [8] [39] . The present work considers the 
following simplified model system of a radiating gas 

Ut + f{u) x + Lq x = 0, 
-q xx + q + M(u) x = 0, 

consisting of a single regularized conservation law coupled with a scalar elliptic equa- 
tion. In (|l.ip . (x, () Six [0, +oo), u and q are scalar functions of (x, t), L G R is a 
constant, and /, M are scalar functions of u. Typically, u and q represent velocity and 
heat flux of the gas, respectively. When the velocity flux is the Burgers flux function, 
f{u) = \u 2 \ and the coupling term M(u) = Mu is linear (M constant), this system 
constitutes a good approximation of the physical Euler system with radiation [5] , and 
it has been extensively studied by Kawashima and Nishibata [THJ [17l [18] , Serre [37] 
and Ito [13] . among others. For the details of such approximation the reader may 
refer to PI1Q21E]. 

Formally, one may express q in terms of u as q — —ICM(u) x , where /C = (1 — 
d x )~ x , so that system represents some regularization of the hyperbolic (inviscid) 
associated conservation law for u. Thus, a fundamental assumption in the study of 
such systems is that 

L^£(u)>0, (1.2) 
au 
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for all u under consideration, conveying the right sign in the diffusion coming from 
Chapman-Enskog expansion (see [31)]). 

We are interested in traveling wave solutions to system (jl.lj) of the form 

{u,q)(x,t) = {U,Q){x-st), (U, Q)(±oo) = (tt±, 0), (1.3) 

where the triple (u+,u_,s) is a shock front of Lax type of the underlying scalar 
conservation law for the velocity, 

«t + /(«)» = 0, (1.4) 

satisfying Rankinc-Hugoniot condition f{u + ) — f(u—) = s(u + — u_), and Lax entropy 
condition ^(u+) < s < 4£(it_). Morover, we assume genuine nonlinearity of the 
conservation law (|1 .4[) . namely, that the velocity flux is strictly convex, 



l_f 
du 



for all u under consideration, for which the entropy condition reduces to u+ < U-. 
We refer to weak solutions of the form (|1.3[) to the system (jl.ip , under the Lax shock 
assumption for the scalar conservation law, as radiative shock profiles. The existence 
and regularity of traveling waves of this type under hypotheses (| 1 . 2|) is known [TBI [22] , 
even for non-convex velocity fluxes [22] . 

According to custom and without loss of generality, we can reduce to the case of a 
stationary profile s = 0, by introducing a convenient change of variable and relabeling 
the flux function / accordingly. Therefore, and after substitution, we consider a 
stationary radiative shock profile (U,Q)(x) solution to (jl.lj) . satisfying 

/([/)' + LQ' = 0, 
-Q" + Q + M(U)' = 0, 1 ' ' 

(here ' denotes differentiation with respect to x), connecting endpoints (u±, 0) at ±oo, 
that is, 

lim (U,Q)(x) = (u±,0). 

X — >±OG 

Therefore, we summarize our main structural assumptions as follows: 



/, M g C 5 
d 2 f 

/(«-) = /(«+) 

r dM I ^ n 
L —r-(u) > 

du 



(regularity), (AO) 

(genuine nonlinearity), (Al) 
(Rankine-Hugoniot condition), (A2) 



(Lax entropy condition), (A3) 
(positive diffusion), (A4) 
where u G For concreteness let us denote 



a(x) := 6(x) := ^(CA(x)), (1.7) 

and assume (up to translation) that a(0) = 0. Besides the previous structural as- 
sumptions we further suppose that 

£6(0) + (k + i)o'(O) > 0, k = 1, . • . , 4. (A5 fe ) 
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Remark 1.1. Under assumption (A4), the radiative shock profile is monotone, 
and, as shown later on, the spectral stability condition holds. Let us stress that, 
within the analysis of the linearized problem and of the nonlinear stability, we only 
need (A4) to hold at the end states u± and at the degenerating value U(0). 

Remark 1.2. Hypotheses (A5&) are a set of additional technical assumptions 
inherited from the present stability analysis (see the establishment of H k energy 
estimates of Section [6] below, and of pointwise reduction bounds in Lemma I3.4[) and 
are not necessarily sharp. It is worth mentioning, however, that assumptions QA5fcD , 
with k — 1,...,4, are satisfied, for instance, for all profiles with small-amplitude 
\u- - u+\, in view of (fl~2f and \U'\ = OQu- - it+|). 

In the present paper, we establish the asymptotic stability of the shock profile 
(U, Q) under small initial perturbation. Nonlinear wave behavior for system (|1.1|) and 
its generalizations has been the subject of thorough research over the last decade. 
The well-posedness theory is the object of study in [21], [T4J, [15l [12] and [2], both for 
the simplified model system and more general cases. The stability of constant states 
[57] , rarefaction waves [TH] IS] , asymptotic profiles [23] 2] [3J for the model system with 
Burgers flux has been addressed in the literature. 

Regarding the asymptotic stability of radiative shock profiles, the problem has 
been previously studied by Kawashima and Nishibata |16j in the particular case of 
Burgers velocity flux and for linear AI — Mu, which is one of the few available 
stability results for scalar radiative shocks in the literature^. In [16 , the authors 
establish asymptotic stability with basically the same rate of decay in L 2 and under 
fairly similar assumptions as we have here. Their method, though, relies on integrated 
coordinates and L 1 contraction property, a technique which may not work for the 
system case (i.e., u € M n , n > 2). In contrast, we provide techniques which may be 
extrapolated to systems, enable us to handle variable ^f(u), and provide a large- 
amplitude theory based on spectral stability assumptions in cases that linearized 
stability is not automatic (e.g., system case, or ^^(u) variable). These technical 
considerations are some of the main motivations for the present analysis. 

The nonlinear asymptotic stability of traveling wave solutions to models in con- 
tinuum mechanics, more specifically, of shock profiles under suitable regularizations 
of hyperbolic systems of conservation laws, has been the subject of intense research 
in recent years (see, e.g., [TU] S3 I2H1 HZl I1H SOI HH HU EH [321 120])- The unifying 
methodological approach of these works consists of refined semigroup techniques and 
the establishment of sharp pointwise bounds on the Green function associated to the 
linearized operator around the wave, under the assumption of spectral stability. A 
key step in the analysis is the construction of the resolvent kernel, together with ap- 
propriate spectral bounds. The pointwise bounds on the Green function follow by the 
inverse Laplace transform (spectral resolution) formula S3] [57J SO] • The main novelty 
in the present case is the extension of the method to a situation in which the eigen- 
value equations are written as a degenerate first order ODE system (see discussion in 
Section [1 .31 below) . Such extension, we hope, may serve as a blueprint to treat other 
model systems for which the resolvent equation becomes singular. This feature is also 
one of the main technical contributions of the present analysis. 



1 The other scalar result is the partial analysis of Serre |38| for the exact Rosenau model; in the 
case of systems, we mention the stability result of 1251 for the full Euler radiating system under zero- 
mass perturbations, based on an adaptation of the classical energy method of Goodman-Matsumura- 
Nishihara [711301. 
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1.1. Main results. In the spirit of EH1 EH EH] , we first consider solutions 
to (jl.ip of the form (u + U, q + Q) , being now u and q perturbations, and study the 
linearized equations of (jl.lj) about the profile (U, Q), which read, 



u t + (a(x)u) x + Lq x = 0, 
-q xx + q + {b(x)u) x = 0, 



(1.8) 



with initial data u(0) = uq (functions a,b are defined in (|1.7|l ). Hence, the Laplace 
transform applied to system (|1.8p gives 



Am + (a(x)u)' + Lq' = S, 
-q" + q+ (b(x)u)' = 0, 



(1.9) 



where source S is the initial data uq. 

As it is customary in related nonlinear wave stability analyses [TJ EH [43J El ESI 
[27l l40j [42] , we start by studying the underlying spectral problem, namely, the homo- 
geneous version of system (|1.9[) : 

(a(x)u)' = —\u — Lq' 

(1.10) 

q" = q+(b(x)u)'. [ ' 

An evident necessary condition for orbital stability is the absence of L 2 solutions to 
(jl.lOP for values of A in {Re A > 0}\{0}, being A = the eigenvalue associated to 
translation invariance. This strong spectral stability can be expressed in terms of the 
Evans function, an analytic function playing a role for differential operators analogous 
to that played by the characteristic polynomial for finite-dimensional operators (see 
[U [351 EJ [43J EH EH SU HOJ [42] and the references therein) . The main property of the 
Evans function is that, on the resolvent set of a certain operator C, its zeroes coincide 
in both location and multiplicity with the eigenvalues of L. 

In the present case and due to the degenerate nature of system (ll.lOp (observe 
that a(x) vanishes at x = 0) the number of decaying modes at ±oo, spanning possible 
eigenfunctions, depends on the region of space around the singularity (see Section [3] 
below, Remark 13. 1|) . Therefore, we define the following stability criterion, where the 
analytic functions D±(X) (see their definition in (|3.32[) below) denote the two Evans 
functions associated with the linearized operator about the profile in regions x ^ 0, 
correspondingly, analytic functions whose zeroes away from the essential spectrum 
agree in location and multiplicity with the eigenvalues of the linearized operator or 
solutions of (fTTTTJt) : 

There exist no zeroes of D±(-) in the non-stable half plane {Re A > 0} \ {0}. (D) 
Our main result is then as follows. 

Theorem 1.3. Assuming (|A0p -( TA5fc[ ), and the spectral stability condition (fDjl . 
then the Lax radiative shock profile (U, Q) is asymptotically orbitally stable. More 
precisely, the solution (u, q) of (|l.ip with initial data uq satisfies 

\u(x, t) - U(x - a(t))\ LP < C(l + i)-* (1 ~ 1/p) l"oUirw4 
\u(x,t) - U(x - a{t))\ m < C(l + ir 1/4 K| L inH4 

and 

\q(x,t) - Q(x - a(t))\ w i, P < C(l +t)-^ 1 - 1/p) \u a \ L i nH , 
\q(x,t) - Q{x - a(t))\ H * < C(l + i) _1/4 |^oki n ^ 
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for initial perturbation u := t2 — U that are sufficiently small in L 1 R H A , for all 
p> 2, for some a(t) satisfying a(0) = and 

\a(t)\ < C\uo\ L i nH A, \a(t)\ < C(l + t)- l ' 2 \u \ LW , 

where ' denotes the derivative with respect to t. 

Remark 1.4. The time-decay rate of q is not optimal. In fact, it can be improved 
as we observe that |<?(£)|l2 < C\u x (t)\ L i and lu^t)!^ is expected to decay like i" 1 / 2 ; 
we omit, however, the details of the proof. 

The second result of this paper is the verification of the spectral stability condition 
(|Dj) under particular circumstances. 

Proposition 1.5. The spectral stability condition (|Dj) holds under either 

(i) b is a constant; or, 

(ii) \u+ — u—\ is sufficiently small. 
Proof. See Appendix IBID 

1.2. Discussions. Combining Theorem 11.31 and Proposition 11.51 we partially 
recover the results of [16] for the Burgers flux and constant L, M, and at the same 
time extend them to general convex flux and quasilinear M. We note that the stability 
result of [16] was for all smooth shock profiles, for which the boundary (see [16) . Thm. 
1.25(ii)(a)) is the condition LM = 1 = — a'(0); that is, their results hold whenever 
LM + a' (0) > 0. By comparison, our results hold on the smaller set of waves for which 
LM + (9/2)a'(0) > 0; see Remark 11.21 By estimating high-frequency contributions 
explicitly, rather than by the simple energy estimates used here, we could at the 
expense of further effort reduce these conditions to the single condition 

LM + 2a'(0) > (1.11) 

used to prove Lemma T3.4I Elsewhere in the analysis, we need only LM + a'(0) > 0; 
however, at the moment we do not see how to remove (jl.llj) to recover the full result 
of [16] in the special case considered there. The interest of our technique, rather, is 
in its generality — particularly the possibility to extend to the system case — and in 
the additional information afforded by the pointwise description of behavior, which 
seems interesting in its own right. 

1.3. Abstract framework. Before beginning the analysis, we orient ourselves 
with a few simple observations framing the problem in a more standard way. Consider 
now the inhomogeneous version 

u t + (a(x)u) x +Lq x =tf, 
-q X x + q+{b{x)u) x =ip, 

of (|1.8p . with initial data u(x,0) = uq. Defining the compact operator JC := {—d x + 
1) _1 and the bounded operator 

Ju:= -Ld x Kd x (b(x)u), 

we may rewrite this as a nonlocal equation 



u t + (a(x)u) x + Ju = ip - Ld x (JCip) , 
u(x, 0) = uq(x) 



(1.13) 



G 
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in u alone. The generator C :— —{a[x)u) x — Ju of (|1.13p is a zero-order perturbation 
of the generator — a(x)u x of a hyperbolic equation, so it generates a C° semigroup 
e ct and an associated Green distribution G(x,t;y) :— e ct S y (x). Moreover, e ct and G 
may be expressed through the inverse Laplace transform formulae 



1 />7+ioo 
Ct _ 1 / At 



e 

2m 



/ e At (A - C)- 1 d\ 

J 7 — zoo 



1 7 — zoo 

G(x,t;y) = ^- I e xt G x (x,y)d\, 



1 



(1.14) 



'7 — zoo 

for all 7 > 70 (for some 7 > 0), where G\(x,y) := (A — C)~ 1 S y {x) is the resolvent 
kernel of C. 

Collecting information, we may write the solution of (|1.12p using Duhamcl's prin- 
ciple/variation of constants as 

/+00 
G(x,t;y)u (y)dy 
-OO 

ft /* + OC 

G(x,t- s;y)(ip - Ld x (JC ip))(y, s) dyds, 

lO J~oo 

q(x,t) = - d x (b(x)u))(x,t), 

where G is determined through (|1.14|) . 

That is, the solution of the linearized problem reduces to finding the Green kernel 
for the w-equation alone, which in turn amounts to solving the resolvent equation for 
C with delta-function data, or, equivalently, solving the differential equation (|1.9[1 
with source S = 5 y {x). We shall do this in standard fashion by writing ()1.9(1 as 
a first-order system and solving appropriate jump conditions at y obtained by the 
requirement that G\ be a distributional solution of the resolvent equations. 

This procedure is greatly complicated by the circumstance that the resulting 3x3 
first-order system, given by 



(@(x)W) x =A(x,\)W where Q(x) 



a(x) 
I 2 



is singular at the special point where a(x) vanishes. However, in the end we find 
as usual that G\ is uniquely determined by these criteria, not only for the values 
Re A > 70 > guaranteed by C°-semigroup theory/energy estimates, but, as in the 
usual nonsingular case 0, on the set of consistent splitting for the first-order system, 
which includes all of {Re A > 0} \ {0}. This has the implication that the essential 
spectrum of £ is confined to {Re A < 0} U {0}. 

Remark 1.6. The fact (obtained by energy-based resolvent estimates) that £ — A 
is coercive for Re A > 70 shows by elliptic theory that the resolvent is well-defined 
and unique in class of distributions for Re A large, and thus the resolvent kernel may 
be determined by the usual construction using appropriate jump conditions. That 
is, from standard considerations, we see that the construction must work, despite the 
apparent wrong dimensions of decaying manifolds (which happen for any Re A > 0). 

To deal with the singularity of the first-order system is the most delicate and 
novel part of the present analysis. It is our hope that the methods we use here may 
be of use also in other situations where the resolvent equation becomes singular, for 
example in the closely related situation of relaxation systems discussed in [2E1 • 
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Plan of the paper. This work is structured as follows. Section [2] collects some 
of the properties of radiative profiles, and contains a technical result which allows us 
to rigorously define the resolvent kernel near the singularity. The central Section [3] 
is devoted to the construction of the resolvent kernel, based on the analysis of solu- 
tions to the eigenvalue equations both near and away from the singularity. Section 
0] establishes the crucial low frequency bounds for the resolvent kernel. The follow- 
ing Section [5] contains the desired pointwise bounds for the "low-frequency" Green 
function, based on the spectral resolution formulae. Section [S] establishes an auxiliary 
nonlinear damping energy estimate. Section [7] deals with the high-frequency region 
by establishing energy estimates on the solution operator directly. The final Section 
[H] blends all previous estimations into the proof of the main nonlinear stability re- 
sult (Theorem (|1.3|l ). We also include three Appendices, which contain, a pointwise 
extension of the Tracking lemma, the proof of spectral stability under linear cou- 
pling or small-amplitude assumptions, and the monotonicity of general scalar profiles, 
respectively. 

2. Preliminaries. 

2.1. Structure of profiles. Under definition (|f ,7p . we may assume (thanks to 
translation invariance; see Remark lC.5l below) that a(x) vanishes exactly at one single 
point which we take as x = 0. Likewise, we know that the velocity profile is monotone 
decreasing (see 22, 23, 38 or Lemma [C.4I below), that is U'(x) < 0, which implies, 
in view of genuine nonlinearity (|1.5|) . that 

a'(x)<0 Viel and ia(a;)<0 Vie ^ 0. 

From the profile equations we obtain, after integration, that 

LQ = f(u ± )-f(U) >0, 

for all x, due to Lax condition. Therefore, substitution of the profile equations (|f .6|) 
yields the relation 

(a'(x) + L b(x)) U' = -LQ - a(x)U", 

which, evaluating at x — and from monotonicity of the profile, implies that 

a'(0) + L6(0) > 0. (2.1) 

Therefore, the last condition is a consequence of the existence result (see Theorem 
IC.3l below). and it will be used throughout. Notice that (A5i) implies condition (|2.ip . 

Next, we show that the waves decay exponentially to their end states, a crucial 
fact in the forthcoming analysis. 

Lemma 2.1. Assuming (|A0|) - (|A4|). a radiative shock profile (U,Q) of (|l.ip 
satisfies 

(d/dx) k (U -u±,Q) <<7e-"l s l, fe = 0,...,4, (2.2) 

as \x\ — ^ +oo, for some r\ > 0. 

Proof. As |x| — > +oo, defining a± — a(±oo) and b± = &(±oo), we consider the 
asymptotic system of (11.61) . that is the constant coefficient linear system 

a±U' = —LQ', 
-Q" + Q = -b±U\ 
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which, by substituting U' into the second equation, becomes 

-Q" -^Q' + Q = 0, 
a± 

or equivalently, 

which then gives the exponential decay estimate (|2.2[) for Q by the hyperbolicity of 
the matrix Aq, that is, eigenvalues of Aq are distinct and nonzero. Estimates for U 
follow immediately from those for Q and the relation 

LQ = f(u±) - f(U), 
obtained by integrating the first equation of (|1.6|) . □ 

2.2. Regularity of solutions near x = 0. In this section we establish some 
analytic properties of the solutions to system (|1.10[) near the singularity, which will 
be used during the construction of the resolvent kernel in the central Section [3] below. 
Introducing the variable p :— b(x)u — q' , system (ll.lOp takes the form of a first-order 
system, which reads 

a(x)u' = - (A + a'(x) + Lb(x)) u + Lp, 

q' = b{x)u-p, (2.3) 
p' = -q. 

For technical reasons which will be clear from the forthcoming analysis, in order to 
define the transmission conditions in the definition of the resolvent kernel (which is 
defined as solutions to the conservative form of system (|2.3j) in distributional sense 
with appropriate jump conditions; see Section T3.ll below), we need p and q to be 
regular across the singularity x = (having finite limits at both sides), u to have (at 
most) an integrable singularity at that point, namely, that u S Lj oc near zero (away 
from zero it is bounded, so this is trivially true), and that it verifies a{x)u — > as 
x — > 0. These properties are proved in the next technical lemma. 

Lemma 2.2. Given A e C, set v := ReA + a'(0) + £6(0)/|a'(0)|. Under assump- 
tions (|A0l) - (| A4[) . and ReA > —Lb(0), then any solution of (|2.3|) verifies 

1. \u{x)\ < C \x\" for x ~ and for some C > 0; 

2. q is absolutely continuous and p is C 1 (for x ~ 0), 

In particular, u € Lj oc (for x ~ 0) and a{x)u{x) ^0 as x — > 0. 

The proof will be done in two steps: (i) first, taking into account "elliptic regu- 
larity" in the equation for p, 

-p"+p = b(x)u, (2.4) 

we prove the Lj oc bound for u close to zero and the subsequent regularity for p and 
q; and (ii), using such a bound, we then prove the pointwise control given in (i). 

Alternatively, one can explicitly solve the above elliptic equation for p and get 
directly the pointwise result for u by plugging the relation into the Duhamcl formula 
for u. Finally, such a control gives the Lj oc property for u and all other regularity 
properties. 
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Proof. [Proof of Lemma |2"T2"] Let us consider the case x > 0, the case x < being 
similar. Consider a fixed x > 0, to be chosen afterwards and let (u,q,p) be any 
solution of (|2.3|) emanating from that point. Therefore, from (|2.4p we know that 



p(x)=Cie X + C 2 e x + / g{x 1 y)u{y)dy 



(2.5) 



for a given (regular) kernel g(x,y). Therefore there exists a constant C Xo such that 
for any e > 

\p\L°°{e, Xo ) < C xo (l + |«|ii( e ,x Q ))< 

Note that the constant C Xg is uniform on e and it stays bounded as xq approaches 
zero and it depends only on the initial values p(xq), q(xo) — —p'(xo). Moreover, the 
Duhamcl principle gives for any x £ [e,£o] : 



u(x) — u(xq) exp 
1 



X + a'(y) + Lb(y) 



dy 



L 



exp 



a (y) 

A + a'(z) + Lb{z) 
a(z) 



dz p{y)dy. 



(2.6) 



From (|1.5p we obtain 

X + a'(x) + Lb(x) X + a'(0) + Lb(0) 



a{x) 

Hence, for x ~ 0, 

r* X + a'(y)+Lb(y) 



a'(0)x 



for x ~ 0. 



exp 



exp 

a; 



A + a'(0) +L6(0) 



A + a'(0) + -C- Kg) 
a'(0) 



a'(0)y 



(2.7) 



Hence the first term of (|2.6[) is integrable in [0,xo] provided Re A > — L6(0), being 
a'(0) < (our argument applies for A ^ -L6(0) - <z'(0); for A = -Lb(0) - a'(0) 
all functions in the integrals above are indeed bounded at zero and the proof of the 
lemma is even simpler). Thus, for a constant C Xo as above, 



\u\me,x ) < U{X )C X0 +C Xo (l+ |u|l1( £ ,i )) X 



1 



■ exp 



x Rc A + a' \z) + Lb(z) , \ 




L \a(y)\ 

Now we use again (|2.7I) to estimate the integral term in (|2.8p as follows: 



dydx. (2. 



1 



Hv)\ 



exp 



Re A + a'(z) + Lb(z) 



xo \a'(0)y\ 



a(z) 
dy = 



dz j dy 
a'(0)\x v 



ReA + a'(0) + Lb{0) 

x 



ReA + a'(0) + Lb{0) 



1 - 



,r 



10 



C.LATTANZIO, C.MASCIA, T.NGUYEN, R.G.PLAZA AND K.ZUMBRUN 



Therefore, 



M.Li(e,xo) < u(x )C Xo + C Xo (l + |u|ii( e , xo ))xo 



»a+H^( W ))»o ,> BeA+L6(0) 



X, 



Re A + L b(0) 







= u(a;o)C Xo + C Xo (l + \u\ L i( €iXO ))x . 



Finally, for a sufficiently small, but fixed xq > 0, from the above relation we conclude 



uniformly in e, namely, u G L 1 (0, eo) for eo > 0. At this point, part 2. of the lemma 
is an easy consequence of expressions (|2.5p , (|2.3|) -j and (|2.3[h . 

Once we have obtained the Lj oc property of u at zero, we know in particular 
\p\l°°(o,x ) is bounded. Hence we can repeat all estimates on the integral terms of 
(|2.6p to obtain part 1 . of the lemma. Finally, 



is again a consequence of Re A > — Lb(0). □ 

Remark 2.3. From condition (|2.1[) it is clear that, for Re A < 0, but sufficiently 
close to zero, u(x) is not blowing up for x — > 0, but it vanishes in that limit, regardless 
of the shock strength (the negative term a'(0) approaches zero as the strength of the 
shock tends to zero). 

3. Construction of the resolvent kernel. 

3.1. Outline. Let us now construct the resolvent kernel for C, or equivalently, 
the solution of (|2.3[) with delta- function source in the u component. The novelty in 
the present case is the extension of this standard method to a situation in which the 
spectral problem can only be written as a degenerate first order ODE. Unlike the 
real viscosity and relaxation cases [26l [27l EH [29] (where the operator C, although 
degenerate, yields a non-degenerate first order ODE in an appropriate reduced space), 
here we deal with the resolvent system for the unknown W := (u,q,p) T 



\u{x)\ L i^ Xo) < a 



lim a(x)u{x) = 



(O(x)W)' = A(a:,A)W, 



(3.1) 



where 




(X + Lb{x)) 
b(x) 






that degenerates at x — 0. 

To construct the resolvent kernel Q\ = Q\(x,y), we solve 

d x (Q(x) g x )-A(x,X) g x = s v {x)i, 



(3.2) 



in the distributional sense, so that 



d x (e(x) g x )-A(x,\)g x = o } 



for all x ^ y with appropriate jump conditions (to be determined) at x = y. The first 
element in the first row of the matrix- valued function g x is the resolvent kernel G x of 
C that we seek. 
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3.2. Asymptotic behavior. First, we study the asymptotic behavior of solu- 
tions to the spectral system 

a{x)u = — (A + a'(x) + L b(x))u + Lp, 

q' = b(x)u-p, (3.3) 

p' = -q, 

away from the singularity at x = 0, and for values of A ^ 0, Re A > 0. We pay special 
attention to the small frequency regime, A ~ 0. Denote the limits of the coefficients 
as 

a± := lim a(x) — ^-(u±), b± := lim b(x) = —r—(u±). 

From the structure of the wave we already have that a+ < < a_ . The asymptotic 
system can be written as 

W'=A±(\)W, (3.4) 

where 

/-a ± 1 (X + Lb±) ag 1 ^ 
A ± (A) := b± -1 

\ -10 

To determine the dimensions of the stable/unstable eigenspaces, let A £ K+, A — > +oo. 
The characteristic polynomial reads 

tt±(a*) := |/u/-A±(A)| =^ 3 + a± 1 (A + J L6±)/i 2 -^-a ± 1 A, 

for which 



3n 2 + 2a± 1 (\ + Lb±)fi-l, 



has one negative and one positive zero, regardless of the sign of a±, for each A > 1; 
they are local extrema of n±. Since n± — > ±oo as /i — > ±oo, 7r±(0) = — a^X has the 
opposite sign with respect to a± and 

7r±(-a±A) =a± (a 2 ± + -^j A 3 + o(A 3 ) A oo, 

so that 7r_/7r + is positive/negative at some negative/positive value of fi, we get two 
positive and one negative zeroes for 7r+ , and two negative and one positive zeroes for 
7r_, whenever A 6 M+, A>1. 

We readily conclude that for each Re A > 0, there exist two unstable eigenvalues 
£ti~(A) and /4(A) with Re/i > 0, and one stable eigenvalue /4(A) with Re/i < 0. 
The stable S + (X) and unstable U + (X) manifolds (solutions which decay, respectively, 
grow at +oo) have, thus, dimensions 

dimf7+(A) = 2, dimS' + (A) = 1, 

in Re A > 0. Likewise, there exist two unstable eigenvalues /uj" (A) , /Uj (A) with Re/i < 
0, and one stable eigenvalue M3 (A) with Re/i > 0, so that the stable (solutions which 
grow at —00) and unstable (solutions which decay at —00) manifolds have dimensions 

dimt/-(A) = l, dim S- (A) = 2. (3.5) 
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Remark 3.1. Notice that, unlike customary situations in the Evans function 
literature [TJ 231 O [23 133 , here the dimensions of the stable (resp. unstable) 
manifolds S + and S~ (resp. U + and U~) do not agree. Under these considerations, 
we look at the dispersion relation 

7T±(if) = -if ~ ait 1 (A + Lb ± )f -i£-a ± 1 \ = 0. 

For each ^£l, the A-roots of last equation define algebraic curves 

A±(0 = -ia±e-^^-, Cei, 

touching the origin at £ = 0. Denote A as the open connected subset of C bounded 
on the left by the two curves A±(£), (eM. Since Lb± > by assumption (A4), the 
set {Re A > 0, A ^ 0} is properly contained in A. By connectedness the dimensions 
of U ± (X) and S ± (X) do not change in A £ A. We define A as the set of (not so) 
consistent splitting 1 , in which the matrices A±(A) remain hyperbolic, with not 
necessarily agreeing dimensions of stable (resp. unstable) manifolds. 

In the low frequency regime A ~ 0, we notice, by taking A = 0, that the eigenvalues 
behave like those of A± (0) . If we define 



9+ := i( - a^Lb+ + Ja+ 2 L 2 b 



= h(aZ X Lb+ + \laZ 2 L 2 b 2 _ 



'l ■— 2 



\{a+ 1 Lb + + Ja- + 2 m\+A) 1 



6»3 := \{ ~aZ 1 Lb++ s JaZ 2 L 2 b 2 _+4), 

as the decay /growth rates for the fast modes (notice that 9j > 0, j = 1, 3), then the 
latter are given by 

/*a (0) = 0, 

fi7(0) = -97 <0<9+ = fi+(0), 
/4(0) = -0+ < < 3 - =Ms(0). 

The associated eigenvectors are given by 

(b ± \l~^{0) 2 f 
V t = ( 

Since the highest order coefficient of 7r± as a polynomial in \i is different from zero, 
then A = is a regular point and whence, by standard algebraic curves theory, there 
exist convergent series in powers of A for the eigenvalues. For low frequency the 
eigenvalues of A± (A) have analytic expansions of the form 

/4(A) = -— + 0(|A| 2 ), 
a± 

nf(\) = ±ef + o(\\\), (3 - 6) 
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corresponding to a slow varying mode and two fast modes, respectively, for low fre- 
quencies. By inspection, the associated eigenvectors can be chosen as 

(b ± \l-^{Xff 
Vf=( -^(A) |. (3.7i 

Notice, in particular, that for this choice of bases, there hold, for A ~ 0, 



V±(X) = [0(A) | , Vf(A) = 0(l), j — 1,3. 



Lemma 3.2. Under the same assumptions as in Theorem \1.3\ for each A € A, 
the spectral system (|3 .4[) associated to the limiting, constant coefficients asymptotic 
behavior of (13 .3|) . has a basis of solutions 




e^W x V±(X), x^O, j = 1,2,3. (3. 
we 

consist of two slow modes 



Moreover, for |A| ~ 0, we can find analytic representations for ^ and , which 



M 2 ± (A) = -a ± 1 A + 0(A 2 ), 

and four fast modes, 

fif(X) = ±ef + o(X), fif(\) = T0f + o(\), 

with associated eigenvectors (|3.7p . 



Proof. The proof is immediate, by directly plugging ()3.8|) into (|3.3p and using the 
previous computations (|3.6[) . (|3.7p . □ 

In view of the structure of the asymptotic systems, we are able to conclude that 
for each initial condition xq > 0, the solutions to (|3.3p in x > xq are spanned by two 
growing modes {ip^(x, A), $2 (a:, A)}, and one decaying mode {$3 (x, A)}, as x — ■> +00, 
whereas for each initial condition xq < 0, the solutions to (|3.3p are spanned in x < xq 
by two growing modes {ip^(x, A), ^(m, A)} and one decaying mode {4>^(x, A)} as 
a; — » — 00. 

We rely on the conjugation lemma of [31] to link such modes to those of the 
limiting constant coefficient system Q3.4p . 



Lemma 3.3. Under the same assumptions as in Theorem ] 1.31 for |A| sufficiently 
small, there exist growing ip^(x,X), j — 1,2, and decaying solutions (f>f(x,X), in 
x 0, of class C 1 in x and analytic in X, satisfying 

1>f(x,\) = e^^V±(X)(I + 0(e-^)), j = 1,2, 
4{x, A) = e^Wy 3 ±(A)(7 + <D(e-*W)), 

where rj > is the decay rate of the traveling wave, and /i^ and Vp are as in Lemma 
HOI above. 

Proof. This a direct application of the conjugation lemma of [3T] (see also the 
related gap lemma in jrJ] in? ] |2"5 ] |2"7] ) . □ 
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As a corollary, and in order to sum up the observations of this section, we make 
note that for A ~ 0, the solutions to (|3.3|) in x > x > are spanned by 

A) = e (^+e(W))* V+(\)(I + 0(e-"l x l)), (fast growing), (3.9) 

i/4(x, A) = e (-V«++0(W 2 ))* V+(X)(I + 0(e-*>W)), (slowly growing), (3.10) 

4>+{x, A) = e (- e 3 + +o(|A|))* 7 3 +(A)(/ + ©(e-'^l)), (fast decaying). (3.11) 

Likewise, all the solutions for x < xo < comprise the modes 

C(x,A) = e <- fl r+ <W))*vr (A)(/ + 0(e"' )|:c| )), (fast growing), (3.12) 

^2~(x,A) = e ( - A/a - +0(|A|2));r y 2 -(A)(/ + C'(e- , 'l :l; l)), (slowly growing), (3.13) 

03 (a, A) = e ( ^+°(l A D) a: y3-(A)(/ + C'(e-' )|a;| )), (fast decaying). (3.14) 

The analytic coefficients V 7 ± (A) are given by (|3.7[) . 

3.3. Solutions near x ~ 0. Our goal now is to analyze system (13 . 3[) close to 
the singularity x = 0. For concreteness, let us restrict the analysis to the case x > 0. 
We introduce a "stretched" variable £ as follows: fix e > and let 




so that £( £ o) = 0, and £ — > +oo as x — > + . Under this change of variables we get 

t du 1 du 1 

da; a(ir) g?£ a(:r) ' 

after denoting ' = In the stretched variables, system (|2.3|) becomes 

/-w L \ 
W = k(£,\)W where A(£,A):= 56 -6), 

V -a / 

and functions oj, a, & are defined by 

W (0:-A + a '(.T(£))+i6W0), 5(0 := a(x(0), b($ :=b(x@). 

Note that from (|2.1[) . for small frequencies A ~ 0, and choosing < eo -C 1 sufficiently 
small we have the uniform bound 

Reu;(£) - Rew(0) = v-= Re A + a'(0) + L 6(0) > 0, 

for all £ G [0, +oo). In addition, we have 

^ = fi(0(a"(z(0) + = O(|2(0l). 

Next, we apply the transformation Z := LIU where 

/l -L/w\ /l L/uA 

L := 1 and R := L" 1 = 1 . 

\0 1 / \0 1 / 
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Since 

/0 Llu ( /lu 2 \ 

|LR| = |LK| = OQa\), and L = =50(1), 

\0 / 

we obtain a block-diagonalized system at leading order of the form 

Z + aQ{£)Z, (3.15) 



-oj 




where 



/0 L/lu L(a" + Lb')/uj 2 ^ 
6= 6 -1 + Lb/w 
\0 -1 

is uniformly bounded. The blocks —ujI and are clearly spectrally separated and 
the error is of order 0(|a(£)|) -> as £ -> +oo. System ([3TT5]) has the form (|Ali|) 
of Appendix [2] (block-diagonal at leading order) and satisfies the hypotheses of the 
pointwise reduction lemma (see Proposition IA.ll below). In our case, there is no 
dependence on a parameter e, Mi = —loI, Mi = and the pointwise error is <5(£) = 
a(£), with constant spectral gap rj. 

Hence, there exist analytic transformations <f>j(£,\), j = 1, 2, satisfying the point- 
wise bound (|A.4[) . for which the graphs {{Z\, $2(^1))}) {($1(^2), Z 2 )} are invariant 
under the flow of (|3.15[) . We now take a closer look at the pointwise error (|A.4|> . 

Lemma 3.4. For the stretched system and for low frequency A ~ 0, there holds 
|$i(£,A)|<Co(0, 3 = 1,2, (3.16) 

provided that 

£6(0) + 2a'(0) > 0. (3.17) 

Proof. From Proposition lA.il there holds the pointwise bound (|A.4|) . namely 
\^,\)\<C [ e-^-^a{y)dy, 







which in terms of the original variables looks like 

\$j(x,\)\ :=\$Mx),\)\<C [*exp (r, J* ^L) dx. 



Since for z small, a(z) ~ a'(0)z, we get 



£ /> (.£ ^) * = (. - 9 W"<™) 



Co'(0)i Co(i) 



?7 + a'(0) t? + o'(0)- 
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in view of < x < e , and as long as a'(0) + r) > 0. Since rj = Re A + a'(0) + Lb(Q), 
condition (|3 . 1 7f> implies (|3.16|) for small A. □ 

Remark 3.5. Notice that (|3.17|) is a stronger condition than (|2.ip , which is 
inherited by the existence result of [22] or Theorem IC. 31 Notably, this new condition 
(|3.17|) holds if we assume (A52). 

In view of the pointwise error bound (|3. 16|) of order 0(a) and by the pointwise 
reduction lemma (see Proposition lA.il and Remark IA.2I below), we can separate the 
flow into slow and fast coordinates. Indeed, after proper transformations we separate 
the flows on the reduced manifolds of form 

Zi = -wZi + ofavzi, 

(3.18) 

Z 2 = (D(~a)Z 2 . 

Observe that the Z\ modes decay to zero as £ — > +oo, in view of 



/|«W<i2 < -(RcA+±7,)£ 



h K -» o, 



as £ — » +oo. These fast decaying modes correspond to fast decaying to zero solutions 
when x — > + in the original u-variable. The Z 2 modes comprise slow dynamics of 
the flow as x — > + . 

PROPOSITION 3.6. Under assumptions (|A0[) - (IA4[) . and (A5 2 ), i/iere exists 
< eo *C 1 sufficiently small, such that, in the small frequency regime A ~ 0, </ie 
solutions to the spectral system (|3.3[) 2?i (~eoiO) U (0,eo) are spanned by fast modes 



(l + 0(o(i))), ±e ^x^0, 



decaying to zero as x — > 0^, anrf slowly varying modes 
±e ^ a; ^ 0, j = 1,3, 
wzi/i bounded limits as x — > 0* . Moreover, the fast modes defined above decay as 





~O(|fa(,))^0, (3.19) 



osn ± , w/iere z/ := (Re A + a'(0) + L6(0))/|a'(0)|. 

Proof. This is a direct consequence of applying our pointwise tracking lemma 
(Lemma l3.4[) to the reduced system ()3.18p . The claimed estimate ()3.19j) for u follows 
in the same way as done in Lemma 12.21 □ 



3.4. Decaying modes. We next derive explicit representation formulae for the 
resolvent kernel Gx{x, y) using the classical construction in terms of decaying solutions 
of the homogeneous spectral problem, matched across the singularity by appropriate 
jump conditions at x = y. The novelty of our approach circumvents the inconsis- 
tency between the number of decaying modes at ±oo. In this section we describe 
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how to construct all decaying solutions at each side of the singularity with matching 
dimensions. 

Choose eo > small enough so that the representations of the solutions of Propo- 
sition [3]6] hold. We are going to construct two decaying modes , j = 1,2 at +oo, 
and one decaying mode W% at — oo. For that purpose, we choose the decaying mode 
at —oo as 

W-( X ,\):=1^ {X > X) > w , (3.20) 

[(71^1 +73^3 + 72W 2 )(x, A), -e < x < 0. 

where the coefficients jj = 7j(A) are analytic in A and such that W% is of class C 1 
in all x < 0. 

To select the decaying modes at +00, consider 

{0, x > 0, 

w;(x,\), -e Q <x<0, (3.21) 

(KlV^r + K 2^2 + K 3<t>3){x, A), X < -C . 

where is the vanishing at x = solution in (|3 . 6|) (the solution is, thus, continuous 
at x = 0), and the coefficients Kj = Kj(A) are analytic in A, and such that the matching 
is of class C 1 a.e. in x. 



Finally, we define 



'<f>£(x,\), X > e , 



W + (x X) ■= t ( aiZ ^ a3Z 3 f + a^it^)^, A), 0<x<e , 
' ' ' * (0!Zi +03Z3 +fow^)(x,\), ~e o <x<0 

Xh^i +hi>2 +^3<f>3)(x,X), x < -e . 

as the other decaying mode at +00, with analytic coefficients a>j, (3j, Sj in A, and 
is of class C 1 a.e. in x. 

Remark 3.7. A similar definition of two decaying modes W% , at —00 and 
one decaying mode at +00, on the positive side of the singularity, is clearly 
available. See Figure l3~Tl 



3.5. Two Evans functions. We first define two related Evans functions 

D ± (y, A) := det(W+ W? Wf)(y, A), for y % 0, 

where Wf = (uf,qf,pf) T is defined as above (see (|332|) . (j3~2T|) . and (|3~20|) ). 
We first observe the following simple properties of D±. 

Lemma 3.8. For A sufficiently small, we have 

D±(y, A) = -aivr'm det K §) + 0(\\\ 2 ), (3.23) 

V^l ^2 / | A=0 

where [u] = u+ — U- . 

Proof. Let us consider (|3 . 23[) for D_. By our choice, at A = 0, we can take 

Wf(x,Q) = Ws{x,0) = W'(x) (3.24) 
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DJy, D 

y<0 




V 


w* ___________ 






* 









DJy, X) 








y>0 


w; 










w 2 











X 



Fig. 3.1. Two Evans functions: Representation of the decaying modes at ±oo for y ^ 0. The 
first picture (left) considers the case y < 0. The sole decaying mode W-J~ at — oo is represented as the 
fast decaying solution on the left; the decaying modes at +oo are two: the exponentially fast decaying 
solution , and the constructed mode , which is identically zero for x > and matched across 
the singularity to the solution which decays to zero as x — > — in the region (y, 0). This provides a 
full set of decaying modes for y < 0. A symmetric construction for the y > case is depicted in the 
second picture (right). 



where W is the shock profile. By Leibnitz' rule, we first compute 
d x D^{ Vl 0) = det (d x W+, W+, W 3 ~) 

^ ' Ia=o 

+ dct (wf,d x W£,W 3 ) 

* ' U=o 



det [W+,W^,d x W : 



where, by using (|3.24p . the second term on the right hand side vanishes and the first 
and third terms can be grouped together, yielding 



d x D_(y, 0) = det (W+, W+, d x W 3 - d x W+ 

Since Wp{-,\) satisfies (|3.1j) . d\W^(x,Q) — (d x uf, d x qf, d x p^) satisfies 

<d{d\W+)' = A{x,0)d x W+(x,0) + d x A(x,0)W+ (x,0), 
which directly gives 

(ad x u+y = -L(d x q+y-u'. 
Likewise, d x W 3 (x,0) = (d x u 3 , d x q 3 , d x p 3 ) satisfies 

(ad x u 3 )' = -L(d x q 3 )' - u'. 



(3.25) 



(3.26) 
(3.27) 



Integrating equations (|3.26p and (|3.27p from +oo and — oo, respectively, with use of 
boundary conditions d x W^(+oo,0) = d\W$ (— oo, 0) = 0, we obtain 



ad x ui = —Ld x ql — u + u+, 



and 



a d 



\u 3 



= -L8 



a<7 3 



Thus 



a(d x u 3 - d x uf) = -L(d x q 3 - d x qf) - [u]. 



(3.28) 



Meanwhile, since Wj, j = 1,2 satisfy the equation (|3.ip and thus (an)' = —Lq' with 
Wf (+oo, A) = 0, we integrate the latter equation, yielding 



-Lqf 



for j = 1,2. 



(3.29) 
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Using estimates (|3. 291) and (|3.28p , we can now compute the A-derivative (|3.25p of D± 
at A = as 



dxD-(y,0) =dct 



= det 



u+ 


«2 


d\Uz 


- d\u{ 


it 


4 


d\<h 


- dxqf 


pi 






- d\pt 




u i 


d\Uz 


- dxu^ 










[u)/L 




p 2 + 


9\P 3 


- dxpt 



(3.30) 



L _1 [m] det 



Pi Pi 



Applying again relation (|3.28[) . we obtain (|3.23p . 
Similarly, for D+ we obtain 

d x D + {y^) = -L- l [u]Aet( U t f) (3.31) 

\Pl P2 J 

from which the conclusion follows. □ 

Since there are two different Evans functions for y ^ 0, we need to be sure if one 
vanishes to order one (part of the stability criterion), then the other does too. Such 
property, content of the following Lemma, guarantees that pole terms are the same 
on y < and y > 0. 

Lemma 3.9. Defining the Evans functions 

D±{\) :=D±(±1,A), (3.32) 

we then have D + (X) = mD-(\) + 0(|A| 2 ) where m is some nonzero factor. 

Proof. Since W^{x) = W is a nonvanishing, bounded solution of the ODE (|3.1[) . 
we must have W r 1 + (1) = miW} (— 1) for some m\ nonzero. Meanwhile, Proposition 
13.61 gives 



uf\ (\x\ 



0(\x\ v a(x)), 



as x — > 0, where v = (a'(0) + i6(0))/|a'(0)|. Thus, smoothness of a near zero guar- 
antees an existence of t\ , ti near zero such that 



u 2 \ =l u 2_ 

± 



This together with the fact that W 2 are solutions of the ODE (|3.1|) yields 



(An 

Pi) x =-i \P2 



x = l 



for some mi nonzero. Putting these estimates into (|3. 30|) and (|3.31|) and using ana- 
lyticity of D± in A near zero, we easily obtain the conclusion. □ 
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4. Resolvent kernel bounds in low frequency regions. In this section, 
we shall derive pointwise bounds on the resolvent kernel G\{x,y) in low-frequency 
regimes, that is, |A| — > 0. For definiteness, throughout this section, we consider only 
the case y < 0. The case y > is completely analogous by symmetry. 

We solve (|3.2p with the jump conditions at x — y: 



[Gx{,y)\ 



'a(y)- 1 N 
1 
1, 



(4.1) 



Meanwhile, we can write Q\(x,y) in terms of decaying solutions at ±oo as follows 



W+(x, X)C+{y, A) + W+(x, X)C+(y, A), x > y, 
-W^(x,X)C^{y,X), x<y 



(4.2) 



where = {C^ k )k = i^2,3 ar e row vectors. We compute the coefficients Cp by means 
of the transmission conditions (|4.ip at y. Therefore, solving by Cramer's rule the 
system 



{W+ W+ Wf) I c+ 



'a(y)" 1 0^ 
1 
1, 



we readily obtain, 




D_(y,X)- 1 (W+ W+ Wi) ad 



l(l/,A) 



'a(y)" 1 0> 
1 
1, 



where M i denotes the adjugate matrix of a matrix M. For example, 



cUy.x) 


= a(y)~ 


- l D-{y,X)- 1 


4 

pi 


% 
Ps 


(y,x), 




= a (yY 




% 

Pa 


it 

pi 


(y,x), 


C^(y,X) 


= a (yY 


-'D-iy^X)- 1 


Qt 
Pt 


pi 


(y,x)- 



(4.3) 
(4.4) 
(4.5) 



Here, note that these are only coefficients that are possibly singular as y near zero 
because of singularity in the first column of the jump-condition matrix (|4.ip . 
We then easily obtain the following. 

Lemma 4.1. For y near zero, we have 

1 



Cj"(y,A)= ^N _1 (l, -L, 0) + 0(1), 
Cr(tf,A) = ~[«]- 1 (l, —L, 0) + O(l), 



(4.6) 



and 



C+(y,X)=a(y)- 1 \yrO(l), 



(4.7) 
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where v is defined as in Proposition [XH and 0(1) is a uniformly bounded function, 
possibly depending on y and A. 

Proof. It suffices to estimate Cj^ when the singularity plays a role. Recalling 
(|3.23[) and (|4.5[) . we can estimate C^(y, A) as 



C 3 "i(2/,A) = -— det 



<?1~ it 



it it 
pi vt 



(.'/•«') = -JL + 0(l), 



where 0(1) is uniformly bounded since a(y)D-(y, A) and normal modes Wp are all 
bounded uniformly in y near zero. This yields the bound for as claimed. The 
bound for C{j follows similarly, noting that Wg - = W± at A = 0. 

For the estimate on , we first observe that by view of (|3.23[) and the estimate 
(pHl?)) on ut, 

D-{y,\)>c\\y\ v , 

for some c > 0. This together with the fact that W% = at A = yields the 
estimate for as claimed. □ 

PROPOSITION 4.2 (Resolvent kernel bounds as \y\ — > 0). Assume (|A0[) - JA5fc[ ). 

For j/ near zero, i/iere /io/d 

g x (x,y)=\- 1 [u]- 1 W'(l, —L, 0) + O(e-^) (4.8) 

/or y < < x, and 

g x (x,y) = \-i[u]- 1 W'(l, -L, O ) + 0(l)(l + T ^— ) (4.9) 
for y < x < 0, and 

gx(x,y)=\- x [u]- x W'(l, —L, 0) + O(e~^\) 
for x < y < 0, /or some ry > 0. Similar bounds can be obtained for the case y > 0. 

Proof. For the case y < < x, using (|4. 6|) and recalling that W*(x) = W' + 
0(A)e" J)|a:| and W£(x) = 0, we have 

Gx(x,y) = W?(x)C+(y) = (w' + 0{\)e-"W)Q^(l, —L, O) + 0(l)), 

yielding (|4.8p . In the second case y < x < 0, from the formula (|4.2p projected on the 
first component, we have 

C+(y, X)u+(x, A) + C+(y, X)u+(x, A) 

where the first term contributes A _1 [u] _1 W A/ + 0(l) as in the first case, and the second 
term is estimated by (|4.7[> and (|3.19|) . 

Finally, we estimate the last case x < y < in a same way as done in the first 
case, noting that y is still near zero and W^{x) = W + 0(X)e~ r, ' x ^. □ 

Next, we derive pointwise bounds of G\(x, y) in regions \y\ — > +oo. Note however 
that the representations (|4.2[) and above estimates fail to be useful in the y — > — oo 
limit, since we actually need precise decay rates in order to get an estimate of form 

\G x (x,y)\<Ce-*\*-*\ 
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which are unavailable from Wj~ in the y — > — oo regime. Thus, we need to express 
the (+)-bases in terms of the growing modes tf>~ at — oo, and the sole decaying mode 
<p 3 where ip~ , (j) 3 are defined as in Lemma 13.31 Expressing such solutions in the basis 
for y < 0, away from zero, there exist analytic coefficients otjk '■— ayfc(A) such that 

Wf(x,X) = a n (X)ipi (x, X) +a 12 (X)ip^(x,X) + a 13 (X)(f> 3 (x, A), 
Wl-(x,X)=a2i(X)iPi(x,X)+a22(X)ip2(x,X)+a 2 3{X)<f>3(x,X). N ' ' ' 

At A = 0, we choose Wf(-,Q) = (j>j(;0) = W'. Thus, 

an(0)=aia(0) = 0. (4.11) 

Furthermore, still as A = 0, ip2 i s a (nearly constant) bounded solution and has the 
form (&_ ,0, 1) T as x near zero. Observe also that W 2 is the solution converging to 
zero in form of \x\ u (l,a(x),a(x)) T as x — > 0~. Thus, we can choose 



a 22 (0) = 0. 



(4.12) 



To express the coefficients Cj~ in terms of the uniformly decaying/growing modes at 
— oo, with a slight abuse of notation we write 

$7 = («7»9T»i'7) T ' 

and define the 2x2 minors 



1,2, 



P 4 P 



-«£(V,A), 



and the analytic minors 
di2(X) 



an ai2 

CK21 Ot22 



d2s(X) 



ai2 ai3 

«22 «23 



<MA) 



an a i3 

«21 a23 



where by (|4.11|) and (|4.12|) . we note 

di 2 (0) = 4,(0) = 0. 
By elementary column operations we notice that 



it %_ 

vi p 3 



= a 2 i0 13 + a 22 rin 



(4.13) 



(4.14) 



P3 Pi 



(4.15) 



\ 1 - ^12 — ^12^i 2 + ^13^13 + d 2 3^ 23 - 
Pi P2 

Lemma 4.3. The minor in (|4.15[) can 6e improved by 

% q t =X(a 11 ^ 1 + a 12 n^) 
P3 Pi 



(4.16) 



(4.17) 
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for some coefficients &ij . 

Proof. The estimate is clear, due to the fact that at A = 0, we can take W±{x, 0) = 
Wf(x,0) = W'(x). □ 

We also have the following crucial cancelation for x > y 
W+C+ + W+C+ 

= a~ x DZ l ({a 21 VL w + a 22 ^){ a n^i + "12^2" + "13^3 ) 
- (anfii3 + 0:12^3) (aiiipi + a 22 ip 2 + "23^3 



(4.18) 



= d 1 D_ 1 (dx2&23*Pl -^12^13^2 - (^13^13 + ^23^23) ( / ) 3 

where VL~- are functions in y and cf>± , (f>Z , ipZ are in x, noting that flZ^ipi an d ^23^2" 
are canceled out. 

We recall here that /i~ , j = 1, 2, 3, are three eigenvalues satisfying 

Hi < -c < 0, iq = -A/a_ + C(A 2 ), \q > c > 0, 
for some Co > 0. 

Lemma 4.4. For y < away from zero, we have 

Cf(y,\) = \- 1 [u]- 1 e-'* v (l, ~L, 0) + O{e-^ y + e-^ y ), (4.19) 
C+(y, A) = 0(e-^ y + e'^ y ) (4.20) 
CZ(y,X) = -\- 1 [u]- 1 e-^ y (l, —L, 0) + O{e-^ y ). (4.21) 

Proof. First, by using ()4.10|) and the estimates in previous sections on normal 
modes ipj , </> J , we readily obtain 

\D_(y,X)\ = \dct(W+,W+,W~)\ = 0(X)e^ y e^ y e^ y 

and 

l^ia I = 0(e^ y e^ y ), |fii 3 | = 0{e^ v e^ y ), jO^sl = C{e^ y e^ y ). 
Using (14. 14|) and noting that 022(0) = 0, we estimate 

C+(y, A) = -a~ l DZ l (^21^3 + "22^23) = \~ l O{ e -^ y ) + 0{e~^ + e~^ y ), 

which gives (|4.19p for A small, by observing that the coefficient in Laurent expansions 
at order A -1 is [u] _1 (l, —L, 0) (see the proof of Lemma 14. II) . 
Next, by (|4.17p . we estimate 

C+(y, A) = a^DZ 1 A(diiOj a + a 12 nz 2 ) = 0(e~^ y + y ). 
Finally, we can estimate 



C 31 (y,X) = a D_ [di 2 n i2 + di 3 fl 13 + d 23 n~ 23 1 

= \- 1 0{e-^ y ) + 0{e-^ y + e-^ v + e~^ v ), 
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in the same way as done for C^, yielding the estimate as claimed; note that the 
constraint (|4.13[) on d(0) shows that only slow-growing mode tp^ appears in the A^ 1 
term. The appearance of the fast-growing term is due to D~ l Qi 2 , a new feature as 
compared to the estimate of C-j^ . □ 

PROPOSITION 4.5 (Resolvent kernel bounds as \y\ — ► +00). Under (|A0[) - QA5fcD ; 
for \y\ large, there hold 



(4.22) 



(4.23) 



(4.24) 



G x (x,y) =X- 1 [u]- 1 e-^ y W'{l, -L, 0) 
+ 0((e-"^ + e""' : V' 1 ) 

for y < < x, and 

G x (x,y) =X- 1 [u]- 1 e-^yW'(l, -L, 0) 

+ 0(e^ {x - y) ) + 0(e^ (x - y) ) + 0(e-^ y e^ x ) 

for y < x < 0, and 

G x (x,y)=-\- 1 [u]- 1 e-^yW'(l, —L, 0) 
+ 0(e-^ y e^ x ) + 0(e^ {x ~ y) ) 

for x < y < 0. Similar bounds can be obtained for the case y > 0. 

Proof. For the first case y < < x, since W 2 = on (0, +00), we have 

g x (x,y) = W+(x)Cf(y,\) = (w'(x) + 0(\e^ x ))c+(y,\) 
= (w'(x) + 0(\e^ x )) 

x (\- l [u)- l e-^ v {l, -L, 0) + O(e-^ y + e -^ v )\ 

giving the estimate (|4.22[) . 

For the third case x < y < 0, we have 

g x (x,y) = -Wz(x)Cs(y,\) = -(w'(x) + 0(Xe^ x ))c^(y, A) 
= (W'(x) + 0(Xe^ x )) 

x (- A- 1 [u]^ 1 e-^y(i ; -L, 0) + 0(e-^ y + e-^ y + e~^ y )y 

proving the estimate (|4.24p . 

Finally, for the second case y < x < 0, we have 

g x (x, y) = W+(x, A) C+(y, A) + W+(x, A) C+(y, A). 

In this case, besides the fact that the A -1 term comes from the expression as 
above, there is a crucial cancelation as computed in (|4.18|) . which proves (|4.23|) . using 
the crucial constraint (|4. 13[) (by our choice of the bases) to eliminate fast modes in 
the A -1 term. □ 
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5. Pointwise bounds and low-frequency estimates. In this section, using 
the previous pointwise bounds (Propositions 14.21 and 14.5ft for the resolvent kernel 
in low-frequency regions, we derive pointwise bounds for the "low-frequency" Green 
function: 

G 7 (z,t;</):=^ / e xt G x {x,y)d\ (5.1) 

where T is any contour near zero, but away from the essential spectrum, and r > is 
a sufficiently small constant such that all previous computations on G\ hold. 

Proposition 5.1. Assuming (|A0j) - flA5fc| ) and defining the effective diffusion 
Lb± (see \27Jl ), the low-frequency Green distribution G 1 {x,t;y) associated with the 
linearized evolution equations may be decomposed as 

G I {x,t;y)=E + G I + R, 

where, for y < 0: 

E{x,t;y) :=U x (x)[u]- 1 e(y,t), 



e(y, t) :— \ errfn [ —. ^ = | — crrfn [ —. = 



|^G J (z,t;y)| < c 1 t-(l/ 3 l+l«l)/ 2 - 1 /2 e -(-y- Q -*) 2 /c 2 t ; 



R(x, t; y) = O( e -n(k-»l+0) + 0{eT*) X {x, y) [l + ^T^/vY 
for some r\, C\, C^> 0, where (3, k = 0, 1 and v = ^"/(q^^ an d 



1 -l<y<x<Q 
otherwise. 



x(x,y) = 

Symmetric bounds hold for y > 0. 

Proof. Having resolvent kernel estimates in Propositions 14.21 and 14.51 we can now 
follow previous analyses of [431 E3 [27] . Indeed, the claimed bound for E precisely 
comes from the term \~ 1 \u\~ 1 e~^ 2 V U X , where = — A/a_ + 0(\ 2 ). Likewise, 
estimates of G 1 are due to bounds in Proposition 14.51 for y away from zero and those 
in Proposition ^. 21 for y near zero but x away from zero. The singularity occurs only 
in the case — l<y<a;<0, as reported in Proposition 14.21 In this case, using the 
estimate (|4.9|) and moving the contour Y in (|5. 1(1 into the stable half-plane {Re A < 0}, 
we have 



XI I 



1+ ' . . )d\ = 0{e- r > t )(l + 



r V a(y)\y\»J 7 V a{y)\y\ v r 

which precisely contributes to the second term in R(x, t; y). The first term is as usual 
the fast decaying term. □ 
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With the above pointwise estimates on the (low-frequency) Green function, we 
have the following from [351 [57] . 

Lemma 5.2 ([251 [57]). Assuming |[A"0]) - QA5 fc D , & satisfies 



+ 00 



d^G I (;t;y)f(y)dy < £7(1 + t)-^-^-^ 2 \f\ L „ (5.2) 

" LP 

for all t > 0, some C > ; for any 1 < g < p. 
We recall the following fact from [41] . 
Lemma 5.3 ([41). The kernel e satisfies 

le^^OlLpJe^-^)!^,^^-^ 1 - 1 /?), 

\e vt (-,t)\ LP <cr^-^-^. 

for all t > 0, some C > 0, for any p > 1. 

Finally, we have the following estimate on J? term. 
Lemma 5.4. Under (|A0j) - J A5fc D , R(x,t;y) satisfies 

R(;t;y)f(y)dy < Ce^\\f\ LP + \f\ L ~), 

LP 

for all t > 0, some C, 77 > 0, /or any 1 < p < 00. 

Proof. The estimate clearly holds for the fast decaying term e~ r, (' x ~ v ' +t ^ in R. 
Whereas, to estimate the second term, first notice that it is only nonzero precisely 
when —1 <j/<x<0or0<a;<j/<l. Thus, for instance, when —1 < x < 0, we 
estimate 



+00 



x(x,y) 



1 



a{y) 

<C\f\ L o 



{ x /vY f(y)dy 
1 



1 



a (y) 



(x/yY f{y)dy 



1 



Hy)\ 



(x/yYdy <C\f\ 



where the last integral is bounded by that fact that a(x) ~ x as |ar| — > 0. From this, 
we easily obtain 



+30 



1 



a(y) 



Lp(-1,0) 



<Ce-^|/|i«, 



which proves the lemma. □ 



Remark 5.5. We note here that the singular term a~ 1 (y)(x/yY appearing in 
(|4.9p and (|5 . 1 [) contributes to the time-exponential decaying term. Note that this 
part agrees with the resolvent kernel for the scalar convected-damped equation ut + 
au x — —Lgu, for which we can find explicitly the Green function as a convected time- 
exponentially decaying delta function similar to terms appearing in the relaxation or 
real viscosity case. 
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6. Nonlinear damping estimate. In this section, we establish an auxiliary 
damping energy estimate. We consider the eventual nonlinear perturbation equations 
for variables (u, q) 

u t + (a(u)u) x + Lq x = d(U x + u x ), 

(6.1) 

-q xx + Q+ {b{u) u) x = 0, 
where a represents the shock location and 

fi(u) := £(C/ + u) = 4L(U) + 0(\u\) = a(x) + 0(\u\), 

b(u) := ^-(U + U ) = ^(U) + 0(\u\) = b(x) + 0(\u\). 

(see subsequent Section [8}. In view of 

a x = a!{x) + 0(\u\ + \u x \), 

then, under assumptions (| A4|) . QA5fc| ), we get that, for all |u|oo and |wa;|oo sufficiently 
small, there holds 

Lb+ (k+ \)a x > 0, (6.2) 



for all k = 1,2,3,4 and all x ~ 0. We are going to profit from (|6.2j) to prove the 
following 

Proposition 6.1. Assume (jAOp - flA5fc| ). So long as 1^1^2,00 and \a\ remain 
sufficiently small, we obtain 



\u\ 2 Hk (t)<Ce-^\u\ 2 Hk (0) + C / e-"(*- s )(| W |i 2 + |d| 2 )( S )d S , V > 0, (6.3) 

Jo 

fork=l,...,4. 

Proof. Let us work for the case a = 0. The general case will be seen as a 
straightforward extension. For our convenience, we denote the <f>— weighted norm as 

k 

1/1*5 : Y^o.yjjyj) : 

for nonnegative functions cj). Now, taking the inner product of q against the second 
equation in (|6.ip and applying the integration by parts, we obtain 

\Qm\h + \<l\h = fa,**) < l\<lx\h + C W\h- 

In fact, we also can easily get for k > 1 

\q\ H H < C\u\ H k-i, (6.4) 

Likewise, taking the inner product of u against the first equation in (|6.ip and inte- 
grating by parts, we get 



\~dj) u \ 2 L* = ~h J a x \u\ 2 dx - (Lq x ,u) 
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which together with (|6.4[) and the Holder inequality gives 



j t Hh<C\u\l 2 . (6.5) 

In order to establish estimates for derivatives, for each k > 1 and <fi > to be 
determined later, we compute 

l±(d k u ,<j>d* u ) = (d k x u u <t>d k x u) = -(Ld k+1 q + d k+1 (au),<pd k u) (6.6) 
where, using the equation for q, we estimate 

-(Ld k+1 q,<pd k u) = -(Ld k (bu) + Ld k -\^d k x u) 

< -{Lb<pd k u,d k u)+e\d k x u\ 2 L2 + C t \\u\ 2 Hk ^ + 

<~l«+ c l u l^-' 

and 

-(^ +1 (au), 09^} = -(a^ +1 w + (fc + l)a*fl£« + L.O.T,<f>d k u) 

= ((i(a^»)x - (fe + l)M)^«,flJ«> - {L.O.T^d k u) 

< {(\(a<P) x -(k + l)a x d>)d k u, d k u) + e\d k u\l* + C e \u\l^ 

By choosing 4> := |a| 2fc+1 , we observe that 

|(6^)» - (fc + 1)M = o 

and thus 

-(d k+1 {au)^d k u) < e\d k u\% , + C £ \u\ 2 H ^ 

for any positive number e. Taking e small enough and putting these above estimates 
together into (|6.6p . we have just obtained 

j t {d k xUl \a\ 2k+1 d k u) < - m (d k u,\a\ 2k+1 d k u) + C\u\ 2 Hk - l7 (6.7) 

for each k > 1 and some small 9\ > 0. 

In addition, by choosing = 1 in (|6.6[) . we obtain 

~(a*«,^ti> = -((L6+ (fc + i)a.)^ U ,5^) + e\d k u\ 2 2 + C,|«|^-i, (6-8) 

for any e > 0. By assumption (|6.2p . there exist 772 sufficiently small and M > 
sufficiently large such that 

M0i|a| 2&+1 + (Lb +(k + i)a x ) > 772 > 0, (6.9) 

for all x S K (by taking M large enough away from zero; for x ~ the bound follows 
from (|6.2p ). Therefore, by adding (|6.8p with A/ times (I6.7[) . using (|6.9p . and taking 
e = 772/2 in (|6.8p . we obtain 

±((l + M\a\ 2k+1 )d k u,d k u}<-^ l^l^+CH^. (6.10) 



STABILITY OF SCALAR RADIATIVE SHOCK PROFILES 



29 



Now, for 5 > 0, let us define 

k 

£{t) :=Y J S i (^+M\a\* k + 1 )d k x u,d k x u). 

i=0 

Observe that £ (t) ~ \u\ 2 Hk . We then use f|6.5|) and (|6. 10[) for k = 1, ...,4 and take <5 
sufficiently small to derive 

j t S(t)<~ m £(t) + C\u\h(t) 

for some 773 > 0, from which (|6.3p follows by the standard Gronwall's inequality. □ 

7. High— frequency estimate. In this section, we estimate the high-frequency 
part of the solution operator e ct (see (|1.14I0 

1 f — 7l+ioo 

S ^ = ^i X mmM> _^e Xt (X-C)-'dX, (7.1) 

for small constants 71,72 > (here \i is the characteristic function of the set /). 

Proposition 7.1 (High-frequency estimate) . Under assumptions (|A0[) - QA5fcD ; 
we obtain 

\dZS 2 (t)(<t>-Ld x (l&lj))\ L2 <Ce-^*(\il>\ H ^ + \<p\H^) k = 0,1, (7-2) 

for some 771 > 0, where K, = {—d 2 , + 1) and L is a con stant (see ([T33]) ). 

Our first step in proving (|7.2p is to estimate the solution of the resolvent system 

Am + (a(a;) u) x + Lq x = tp, 
-q-xx +Q + (b(x) u) x = tp, 

where a(x) — —(U(x)) and b(x) = ——(£/) as before. 
du du 

Proposition 7.2 (High-frequency bounds). Assuming (|A0j) - QA5fcD , for some 
R, C sufficiently large and 7 > sufficiently small, we obtain 

|(A - C)~ l {ip - Ld x (JCi/}))\ H i <c(\ip\ 2 H1 +\i>\ 2 L2 ), 

I (A - C)-\ V - Ld x (lCtl>))\v < uTTTi (l^" 1 + l^ 2 ) > 
for all |A| > R and Re A > —7. 



Proof. A Laplace transformed version of the nonlinear energy estimates (|6.5p and 
(fF7TU|) in Section E] with k = 1 (see 02], pp. 272-273, proof of Proposition 4.7 for 
further details) yields 

(ReA+^)H^<c(Hl 2 + M^ + H| 3 ). (7.3) 

On the other hand, taking the imaginary part of the L? inner product of U against 
Xu = Cu+d x LICh+f and applying the Young's inequality, we also obtain the standard 
estimate 

|ImA|M£ a < \(Cu,u)\ + \{LK^,u x )\ + \(<p,u)\ 
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noting the fact that £ is a bounded operator from H 1 to L 2 and AC is bounded from 
L 2 to H 1 . 

Therefore, taking 7 = 71 /4, we obtain from (|7.3[) and (|7.4[) 



for any Re A > —7. Now take i? sufhciently large such that \u\ 2 L2 on the right hand 
side of the above can be absorbed into the left hand side for |A| > R, thus yielding 

\\\Mh + Wm <c(Hi 2 + M^), 

for some large C > 0, which gives the result as claimed. □ 
Next, we have the following 

Proposition 7.3 (Mid- frequency bounds). Assuming (|A0p - flA5fcD , we obtain 
|(A - £)~VU a <C\ip\ H i for R- 1 < |A| < R and Re A > -7, 

/or any R and C — C(R) sufficiently large and 7 = 'y(R) > sufficiently small. 

Proof. Immediate, by compactness of the set of frequency under consideration 
together with the fact that the resolvent (A — £) _1 is analytic with respect to H l in 
A; see, for instance, [41]. □ 

With Propositions 17.21 and 17.31 in hand, we are now ready to give: 

Proof. [Proof of Proposition 17.1] The proof starts with the following resolvent 
identity, using analyticity on the resolvent set p{C) of the resolvent (A — £) -1 , for all 
if e £>(£), 

(A - £)~V = A _1 (A - £)- x £^ + A~Y 

Using this identity and (|7.1[) . we estimate 

-71+ioo 



1 r • 
s 2 (t) v = — \ 



71 — 200 

^ /- — 71+ioo 

=-S 1 + S 2 , 

where, by Propositions 17.11 and 17.31 we have 

-71+icxD 



X { ^>^X- x {\-L)- x C V d\ 



|5i 



L- 



< c 



|Ar l e RcA t|(A _ £r l^ |i2|dA| 



71 — too 



< Ce" 71 ' ( 

-71 — too 

< Ce~^W\ m 



|A|- 3 / 2 |£^| Hl |dA| 



and 



15: 



2|L 2 



< 



1 

2^ 



— 7l +200 



A~V*dA 



-71 — 'iOO 

<C7e- 7lt M i2 , 



L' 2 



1 

2^ 



-71+ir 

■71— ir 



A _1 e A *dA 



Z+ 
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by direct computations, noting that the integral in A in the first term is identically 
zero. This completes the proof of the bound for the term involving 93 as stated in the 
proposition. The estimate involving ip follows by observing that Ld x K, is bounded 
from H s to H s . Derivative bounds can be obtained similarly. □ 

Remark 7.4. We note that in our treating the high-frequency terms by energy 
estimates (as also done in HOI 131])) we are ignoring the pointwise contribution there, 
which would also be convected time-decaying delta functions. To see these features, 
a simple exercise is to do the Fourier transform of the equations about a constant 
state. 

8. Nonlinear analysis. In this section, we shall prove the main nonlinear sta- 
bility theorem. Following [TTII27] , define the nonlinear perturbation 

Q (x,t) := Q (x + a(t),t) - (^j (*), (8.1) 

where the shock location a(t) is to be determined later. 

Plugging (|8.ip into (|1.1|) . we obtain the perturbation equation 

u t + (a(x) u) x + Lq x = N x {u) x + a(t) (u x + U x ), 
-Qxx + q + (K x ) u )x = N 2 (u) x , 

where Nj(u) = 0(\u\ 2 ) so long as u stays uniformly bounded. 
We decompose the Green function as 

G(x, t; y) = G'ix, t; y) + G u (x, t; y) (8.2) 

where G I (x,t;y) is the low-frequency part. We further define as in Proposition 15. 11 

G\x, t; y) = G z (x, t; y) - E(x, t; y) - R(x, t; y) 

and 

G II {x 1 t ] y) = G II {x,t-y)+R{x 1 t-y). 
Then, we immediately obtain, from Lemmas l5.2U5.4l and PropositionEH the following 



Lemma 8.1. There holds 

" + OC 

d^G I (;t;y)f(y)dy < C(l + t)-^'"-^-^' 2 \f\ M , (8.3) 



y 

for all 1 < q < p, f3 = 0, 1, and 
r+00 



L'' 



G n (x,t;y)f(y)dy <Ce~^\f\ HS , 



LP 



(8.4) 



for all 2 < p < 00. 

Proof. Bound (|8.3|) is precisely the estimate (|5.2p in Lemma IS~2"| recalled here for 
our convenience. Inequality (|8.4j) is a straightforward combination of Lemma l5.4l and 
Proposition [TTj followed by a use of the interpolation inequality between L 2 and L°° 
and an application of the standard Sobolev imbedding. □ 
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We next show that by Duhamel's principle we have: 
Lemma 8.2. There hold the reduced integral representations: 

/+oo ^ 
(G 1 + G n )(x, t; y)u (y)dy 
-oo 

t r+ co ~ / \ 

/ Gt(x, t - s; y) ( N x (u) - L K d y N 2 {u) + du) (y, s) dy ds 

J-oo V ' 

1 r+°° ~ / \ 

I G II (x,t-s;y)[N 1 (u)-LlCd y N 2 (u) + au) (y,s)dyds, 

q(x,t) =(iCd x )(N 2 {u)-bu)(x,t), 

(8.5) 



and 



r+oo 

*(t) = - / e(y,t)u (y)dy 

(8.6) 

/ e y (y,t - s)lNi(u) - LKd y N 2 (u) + au\(y,s) dyds. 

J-oo v ' 

d(f) = - / e t (y,t)u {y)dy 

(8.7) 

I e yt (y,t - s)[Ni(u) - LKd y N 2 {u) + au)(y,s) dyds. 

J-oo V 7 



Proof. By Duhamel's principle and the fact that 

r+oo 

Ctrril 



G(x, t- y)U'(y)dy = e ct U'(x) = U'(x) 



we obtain 



/+oo 
G(x,t;y)u (y)dy 
-oo 

pi r+oo 



J-oo 

a{t) U'. 



G(x, t — s;y) ^Ni(u) — L/C d y N 2 {u) + a uj (y, s) dy ds 



Thus, by defining the instantaneous shock location: 

/+oo 
e(y,t)u (y)dy 
-oo 
ft /> + oo 



J-oo 



e y (y, t - s) (^N ± (u) - L K d y N 2 (u) + auj(y,s) dy 



ds 



and using the Green function decomposition (18. 2j) . we easily obtain the integral rep- 
resentation as claimed in the lemma. □ 

With these preparations, we are now ready to prove the main theorem, following 
the standard stability analysis of [551 HOI SI] : 
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Proof. [Proof of Theorem II. 3| Define 



((t) := sup 

0<s<t,2<p<oo 



*(*)|£*(1 + S )^ {1 - 1M + \a(s)\ + |d(a)|(l + s) 1 ' 2 



We shall prove here that for all t > for which a solution exists with uniformly 
bounded by some fixed, sufficiently small constant, there holds 

C(t)<C(\u Q \ LlnHs +at) 2 )- (8.8) 

This bound together with continuity of £(t) implies that 

C(*) < 2C\u Q \ L i nH s (8.9) 

for t > 0, provided that \u \Li n H s < 1/4C 2 . This would complete the proof of the 
bounds as claimed in the theorem, and thus give the main theorem. 

By standard short-time theory/local well-posedness in H s , and the standard prin- 
ciple of continuation, there exists a solution u G H s on the open time-interval for 
which \u\hs remains bounded, and on this interval is well-defined and continu- 
ous. Now, let [0, T) be the maximal interval on which \u\h s remains strictly bounded 
by some fixed, sufficiently small constant 5 > 0. By Proposition l6.ll and the Sobolev 
embedding inequality |tt|^2,<» < C|u|jja, s > 3, we have 



l«(*)llr- < Ce-^luolir. + ^ J (|«(r)|| a + |d(r)| 2 )dr 

<c(| Uo |^+CW 2 )(i + *r 1/2 



(8.10) 



and so the solution continues so long as £ remains small, with bound (|8.9p . yielding 
existence and the claimed bounds. 

Thus, it remains to prove the claim (|8.8j) . First by representation (|8.5p for u, for 
any 2 < p < oo, we obtain 



\u\ L p(t) < 



+ 00 



(G 1 + G II )(x 1 1; y)u (y)dy 



LP 



— oo 



G l y (x, t - s; y) (Ni(u) -LJC d y KN 2 (u) + d (y, s) dy 



ds 



L" 



G H (x,t 



;y)(N 1 (u)-LJCd y N 2 (u)+au^(y,s)dy 



ds 



= h+h + h, 
where estimates (|8.3[) and (|8.4|) yield 

h = 



LP 



(G 1 + G n )(x, t; y)u (y)dy 

J 

< C(l + t)-* (1 - 1/p) |«oUi + Ce-*\uo\ H , 
<C-(l + f)-* (1 - 1/p) l«oU*n*» 5 
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and, with noting that L fC d v is bounded from L 2 to L 2 



h = 



+00 



G*(x,t - s;y)\Ni(u) - LJCd y N 2 (u) + auj (y, s) dy 

<C I (t - s)-^/ 2 ^/p)-^ 2 (\u\ l ^ + \a\)\u\ L 2( s )ds 
Jo 

< CC(t) 2 [\t - fl )-i(l/2-l/p)-l/2(l + s )-3/4 ds 
JO 

<cat) 2 {i + t)-i (1 ~ 1/p) , 

and, together with (|8. 10[> . s > 4, 



L" 



h = 



G 



11 



(x,t- s;y)(iVi(tt) - LJCd y N 2 {u) +auj {y,s)dy 



ds 



L'< 



<C e-^-^lN^u) - LKd y N 2 (u) + au\ m {s)ds 
Jo 

<C [ e-^- s \\u\ H s + \a\)\u\ H s(s)ds 



< C(\u \ 2 H3 + C(t) 2 ) J e-^-°\l + s )-id S 

<C(\u \ 2 H3 +C(t) 2 )(l + t)-\ 

Thus, we have proved 

\u(t)\ LP (l + t)^ 1 " 1 ^) < C(\uo\ L i nH . + C(t) 2 ). 

Similarly, using representations (|8.6[) and (|8.7p and the estimates in Lemma 15.31 on 
the kernel e(y,t), we can estimate (see, e.g., [251 

\a(t)\(l + t) 1 ' 2 + \a(t)\ < C(\uo\ L i + C(t) 2 )- 



This completes the proof of the claim (|8.8p , and thus the result for u as claimed. 
To prove the result for q, we observe that JCd x is bounded from L p — > W 1,p for all 
1 < p < co, and thus from the representation (|8.5p for g, we estimate 



and 



\q\w^{t) < C(\N 2 (u)\ LP + \u\ LP )(t) 

^C^L^KC^^Hsil + t)-^ 1 - 1 ^ 



(*) < CM*, (t) < CluoUxnff. (1 + t)~ 1/4 > 
which complete the proof of the main theorem. □ 

Appendix A. Pointwise reduction lemma. Let us consider the situation of 
a system of equations of form 



W x = A e (x,X)W, 



(A.l) 



for which the coefficient A £ does not exhibit uniform exponential decay to its asymp- 
totic limits, but instead is slowly varying (uniformly on a e-neighborhood V, being 
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e > a parameter). This case occurs in different contexts for rescaled equations, such 
as p. 151) in the present analysis. 

In this situation, it frequently occurs that not only A e but also certain of its 
invariant eigenspaces are slowly varying with x, i.e., there exist matrices 



for which L e M e (a;) = I and |LR'| = |L'R| < C5 e (x), uniformly in e, where the point- 
wise error bound 5^ — 5 e (x) is small, relative to 



M £ := L e A e R e (a;) = (x) (A.2) 

and ""' as usual denotes d/dx. In this case, making the change of coordinates W e — 
WZ, we may reduce (|A.1[) to the approximately block-diagonal equation 

Z £ ' = M £ Z e + S € Q e Z e , (A.3) 

where M c is as in (|A.2[) , 9 e (a:) is a uniformly bounded matrix, and S e (x) is (relatively) 
small. Assume that such a procedure has been successfully carried out, and, moreover, 
that there exists an approximate uniform spectral gap in numerical range, in the strong 
sense that 

mincr(ReMi) — maxcr(ReM|) > ^(x), for all x, 

with pointwise gap rf(x) > t]q > uniformly bounded in x and in e; here and elsewhere 
ReiV := i(7V + N*) denotes the "real", or symmetric part of an operator N. Then, 
there holds the following pointwise reduction lemma, a refinement of the reduction 
lemma of |27j (see the related "tracking lemma" given in varying degrees of generality 

in Eg G21 E31 io]). 

Proposition A.l. Consider a system (|A.3|) under the gap assumption (fA")) . 
with O e uniformly bounded in e £ V and for all x. If, for all e £ V, sup :reR (i5 e /77 e ) 
is sufficiently small (i.e., the ratio of pointwise gap r) e (x) and pointwise error bound 
^ e (x) is uniformly small), then there exist (unique) linear transformations $f(a;,A) 
and $2(^7^)7 possessing the same regularity with respect to the various parameters e, 
x, A as do coefficients M e and S e (x)Q <i (x), for which the graphs {{Z\, $2(^1))} an d 
{($1(^2), Z2)} are invariant under the flow of (| A.3[) . and satisfying 

sup 1^1 < Csup(5777 e ). 

R R 

Moreover, we have the pointwise bounds 

\*Z(x)\ <C [ X e-Sy^ {z)dz 8^y)dy, (A.4) 



and symmetrically for $f . 

Proof. By a change of independent coordinates, we may arrange that rj e {x) = con- 
stant, whereupon the first assertion reduces to the conclusion of the tracking/reduction 
lemma of [27] . Recall that this conclusion was obtained by seeking $| as the solution 
of a fixed-point equation 

= T& 2 (x) := [ X T^ x 5\y)Q{& 2 ){y)dy. 

J —OO 
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Observe that in the present context we have allowed S e to vary with x, but otherwise 
follow the proof of [57] word for word to obtain the conclusion (see Appendix C of 
[27] . proof of Proposition 3.9). Here, <3($ 2 ) = 0(1 + I^H 2 ) by construction, and 
|^72/-»x| < Q e -ri(x-y) _ Thus, using only the fact that $ 2 | is bounded, we obtain the 
bound (|A.4[) as claimed, in the new coordinates for which rf is constant. Switching 
back to the old coordinates, we have instead l^^l < Ce~^» n ^ dz , yielding the 
result in the general case. □ 

Remark A. 2. From Proposition lA.il we obtain reduced flows 

r z? = m\z\ + 6%e n + ot 2 nWi, 

on the two invariant manifolds described. 
Appendix B. Spectral stability. 

Consider the eigenvalue system (II . 10[) . Integrating the equations we find the 
zero-mass conditions for u and q, 

udx = 0, ^ qdx = o, 

which allows us to recast system (jl.lOj) in terms of the integrated coordinates, which 
we denote, again, as u and q. The result is 

Xu + a(x) u + Lq = 0, 
-q" + q + b(x) u' = 0. V ' ' 

The following proposition is the main result of this section. 

Proposition B.l. Let (u,q) be a bounded solution of (jB.ip . corresponding to 
a complex number A ^ 0. Then Rc A < provided that at least one of the following 
conditions holds 

(i) b is a constant; 

(ii) |w+ — u_| is sufficiently small. 

Proof. In any case, we can assume b > by redefining q by — q if necessary, still 
preserving the condition Lb > 0. Taking the real part of the inner product of the first 
equation against bu and using integration by parts, we obtain 

ReA|6 1/2 u|2 2 = -Re(abu',u) -Re(Lq',gu) 

= Rc (((a b)'u, u) + (Lq, (b u)')) 

= Rc ({(a b)'u, u) + (Lq, q" - q + b' u)\ 

= Rc ((a' bu,u)- (L q' , q') - (Lq, q) + (a b' u, u) + (L b' q, u)^j 
< (a'bu,u)-±\q\ 2 Hl +C({\a\ + \b'\)\b'\u,u), 

d 2 f 

which proves the proposition in the first case, noting a' = — tt(U) U' < (by mono- 

du z 

tonicity of the profile) and b > 6 > 0. For the second case, observe that \a\ + \b'\ is 



STABILITY OF SCALAR RADIATIVE SHOCK PROFILES 



37 



now sufficiently small and \b'\ and \a'\ have the same order of "smallness", that is, 
of order = 0(\u + — u_|). Thus, the last term on the right-hand side of the 

above estimate can be absorbed into the first term, yielding the result for this second 
case as well. □ 

Appendix C. Monotonicity of profiles under nonlinear coupling. 

In this Appendix we show that radiative scalar shock profiles for general nonlinear 
coupling are monotone, a feature which plays a key role in our stability analysis. Al- 
though the existence of profiles for nonlinear coupling is already addressed in [23] , and 
the monotonicity for the linear coupling case is discussed in [35} H2] > f° r completeness 
(and convenience of the reader) we closely review the (scalar) existence proof of [22] 
and extend it to the nonlinear coupling case, a procedure which leads to monotonicity 
in a very simple way. 

The main observation of this section is precisely that, thanks to assumptions (|A0|) 
and (|A4[) . the mapping u i— > LM(u) is a diffeomorphism on its range [23j . which can 
be regarded as the identity along the arguments of the proof leading to the existence 
result of [22]. Since LM is monotone increasing in [u+,u_], setting M± := M(u±), 
there exists an inverse function H : [LM + , LAI-} — > [u + , u_] such that 

y = LM(u) ^ u = H(y), 

for each u G and with derivative 

dH ( dM y 1 
— =[L—{H{y))) > 0. 
ay \ du J 

Consider once again the stationary profile equations (| 1 . 6[) (after appropriate flux 
normalizations), with (U, Q) (±oo) = (u±,0). Integration of the equation for Q leads 
to JmQ = —[M] = M_ — M + . Let us introduce the variable Z as 

Z :=-L f Q(£)d£ + Lb-, 

such that Z' = —LQ and Z — > LM± as x — + +oo. In terms of the new variable Z the 
profile equations are 

z" = f(uy, 

Z' - Z'" = LM(u)'. 

Integrating las equations, and using the asymptotic limits for Z, we arrive at the 
system 



Z' = f(U) - f(u±) 
Z-Z" = LM(u). 



, ^ (ci) 



We can thus rewrite the ODE for Z as 

Z' = F(H{Z - Z'% (C.2) 



where F(u) := f(u) — f(u±). In view of strict convexity of /, the function F is strictly 
decreasing in the interval and strictly increasing in [u*,it_], with F(u±) = 
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and F(u*) = — m < 0. Hence, F is invertible in those intervals with corresponding 
inverse functions h±, and we look at the solutions to two ODEs, namely, 

Z" = Z - LM(h±(Z% 
Z(±oc) = LM±, Z'(±oo) = 0, 

in their corresponding intervals of existence. Observe that the derivatives of the 
functions h± are given by h' ± = 1/ f'(h±(-j), with /'(«) ^ in [u + ,it») U 
Note that /i + : [— m, 0] — > [«+,«*] and : [— m, 0] — > [ti*,u_], and that h + (h—) is 
monotonically decreasing (increasing) on its domain of definition. 

Following [22] closely, we shall exhibit the existence of a Z-profile solution to 
(|C.2[) between the states LM_ > LM + , for which the velocity profile follows by 
U = H(Z — Z") (see (|U.ip ). In the sequel we only indicate the differences with the 
proofs in Section 2 of [22j , and pay particular attention to the monotonicity properties 
of Z, which leads to the monotonicity of U in Lemma [C.4I below. 

The following proposition is an extension of Propositions 2.2 and 2.3 in [55] to 
the variable G' case. 

Proposition C.l ([22]). (i) Denote Z + = Z + (x) the (unique up to translations) 
maximal solution to 

Z" = Z - LM(h+(Z')), 

with conditions Z(+oo) = Lb+ and Z'(+oo) = 0. Then Z+ is monotone increasing, 
Z' + is monotone decreasing, and Z+ is not globally defined, that is, there exists a point 
that we can take without loss of generality as x — (because of translation invariance) 
such that 

Z+(0) - Z'_l(0) = LM{u»), Z' + (0) = -to < 0. 

(ii) Denote Z_ = Z_(x) the (unique up to translations) maximal solution to 

Z" = Z - LM(h-(Z')), 

with conditions Z(—oo) = Lb^ and Z'(-oo) — 0. Then Z_ and Z'_ are monotone 
increasing, and Z_ is not globally defined, that is, there exists a point that we can 
take without loss of generality as x = (because of translation invariance) such that 

Z_(0) - Z'1(0) = LM(u*), Z'_(0) = — to < 0. 



Proof. We focus on part (i) of the Proposition. The second part is analogous. 
Rewrite the equation for Z + as X' = J+(X) with X = (Z. Z') T and 



J+(X) = 

for which 



Z' 

Z - LM(h+{Z')) 



'0 1 
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in view of h + (Q) = u + , and therefore, the starting point (LM+,0) of the trajectory 
is a saddle point. We focus on the stable manifold as we need Z to be decreasing. 
Follow the trajectory that exits from (LM + , 0) in the lower half plane of the phase field 
(Z, Z'). We claim that Z is strictly monotone decreasing and Z 1 is strictly monotone 
increasing. Suppose, by contradiction, that Z attains a local maximum at x$ G R. 
Then Z'(x ) = and > Z"(x Q ) = (Z - LM{h+{Z')) x=XQ = Z(x ) - LM+, which 
is false. Hence, Z is monotone decreasing and Z' < 0. Now, assume that Z' attains 
a local minimum at x = x$. Then the trajectory Z' = tp(Z) in the phase plane must 
attain a local minimum at the same point, yielding ip'(Z) = and <f"(Z) > 0. Thus, 
at x = xo, 

= ip'(Z) = Z"/Z' = (Z- LM(h+(Z')))/Z' 

and 

ip"{Z) = (d/dZ)((Z LM(h + (Z')))/Z') 

= 1/Z' - (dZ'/dZ) ((Z' + LM'(h + (Z'))h' + (Z')Z' - LM(h+(Z')))/(Z') 2 ) . 

But (dZ'/dZ) = ip'(Z) = at x = x , thus <p"{Z) = 1/Z' < 0, which is a contra- 
diction. This shows that Z' is strictly monotone increasing with Z" > 0, and clearly 
LM{h+{Z')) e [LM+,LM{u*)\, h+{Z') e [u +) u*]. This shows that Z" = Z + 0(1) 
and the solution does not blow up in finite time. 

By following the proof of Proposition 2.2 in |22] word by word from this point 
on, it is possible to show that the solution reaches the boundary of definition of the 
differential equation at a finite point which, by translation invariance, we can take as 
x = 0. Hence, Z' + (0) - Z^(0) = LM{u*) and Z' + (0) = -m < hold. This concludes 
the proof. □ 

Lemma C.2. For the maximal solutions Z± of Provosition [CA[ there holds 
Z-(0) < LM{u*) < Z+(0). 

Proof. This follows by mimicking the proof of Lemma 2.4 in |22j . We warn the 
reader to now consider the dynamical system 

y' = F(H{y)), 
y(±oo) = LM±. 

A comparison of the solution y of the system above with the trajectory Z + in the 
phase space yields the inequality on the right. The other inequality is analogous. See 
[2"2] for details. □ 

The last lemma guarantees the existence of a point of intersection for the orbits 
of the maximal solutions Z+ and Z- in the phase state field. The monotonicity of 
Z± and Z± implies that the intersection is unique. Matching the two trajectories at 
that point provides the desired Z-profile. Hence, we have the following extension of 
the existence result in (Theorem 2.5). 

Theorem C.3 (|22j). Under assumptions, there exists a (unique up to trans- 
lations) Z-profile of class C 1 with Z(±oo) = LM±, solution to (|C.2|) . The solution 
Z is of class C 2 away from a single point, where Z' has at most a jump disconti- 
nuity. Moreover, there exists a (unique up to translations) velocity profile U with 
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U(±oo) = u± solution to (|C.1[) . which is continuous away from a single point, where 
it has at most a jump discontinuity satisfying Rankine-Hugoniot conditions and the 
entropy condition. 

Proof. Lemma IC.2I implies the existence of a point in the (Z, Z') plane where 
the graphs of Z_ and Z + intersect. By monotonicity of the graphs the intersection is 
unique. Thus, after an appropriate translation, we can find a point iGl such that 
{Z-(x), Z'_(x)) = (Z + (x), Z' + (x)) =: (Z,Y), and the Z-profile is defined as 

ry, N )Z+(X), X > X, 

Z{ - x > := 7 M S - 
I Z-{x), x < x. 

Z is C 1 and satisfies Z — > Lb± as x — > ±00. Moreover, Z is C 2 except at x = x. The 
velocity profile is now defined via 

U := H(Z-Z"), 

with the described regularity properties due to regularity of Z and the fact that H = 
(iM)" 1 is of class, at least, C 2 . Likewise, at the only possible discontinuity a; = a; of U 
is is possible to prove that U satisfies Rankine-Hugoniot condition, U{x — 0) = U (x+0) 
and the entropy condition U(x - 0) = h^(Y) > h+(Y) = U(x + 0). □ 

Lemma C.4 (Monotonicity). The constructed profile U is strictly monotone de- 
creasing. 

Proof. Let X2 > X\, with Xi ^ x, and suppose that U(x2) > U{x\), that is, 
H(Z — Z")\ X=X2 > H(Z — Z")\ x=Xl . Since H is strictly monotone increasing we 
readily have that 

LM(h±(Z' ± (x 2 )) = (Z- Z'%= X2 >(Z- Z") X=X1 = LM(h ± (Z' ± ( Xl )), 

where the ± sign depends on which side of x = x we are evaluating the Z-profile. Sup- 
pose x\, X2 are on the same side, say, x < x\ < X2 (the symmetric case, x\ < X2 < x, 
is analogous). Since LM is monotonically increasing, last condition implies that 
h+(Z' + {x2)) > h+(Z' + {xi)) . But this is a contradiction with the fact that Z' + is mono- 
tone increasing and h + is strictly decreasing, yielding h+(Z' + (x2)) < h + (Z' + (xi)). The 
case xi < x < X2 leads to the condition h + (Z' + {x2)) > h-(Z'_(xi)), which is obvi- 
ously false in view that h + : [— m,0] — > and h- : [— m, 0] — > yielding 
again a contradiction. Finally, we remark that at the only point of discontinuity of 
[/, namely at x = x, the jump is entropic, satisfying U(x — 0) > U(x + 0). Therefore 
U is strictly monotone decreasing in all 16I. D 

Remark C.5. Observe that the constructed velocity profile is continuous, except, 
at most, at one point where it observes an entropic jump. The regularity of U increases 
as long as the strength of the profile decreases below an explicit threshold [TBI HI] > 
becoming continuous and, moreover, of class C 2 . We remark, however, that away 
from the possible discontinuity x = x, the profile has the same regularity of Z", 
independently of the shock strength, because of smoothness of H. Whence, from 
regularity assumption (jAOjl and by differentiating equation (|C.3j) . Z" is of class C 2 
away from x — x, and so is U . Finally, thanks to translation invariance we have chosen 
x = to be the point where the equations for the profiles Z± reach LM(u*), being 
the only zero of £(it); this implies that U(0) = H(Z-Z") X=0 = H(LM(h±(-m))) = 
u», so that a(x) — ^{U) vanishes only at x = 0. 
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In view of last remark we have the following 

Corollary C.6. Except for a possible single point x = x, the profile U is of 
class C 2 and satisfies U' < a.e. Moreover, the function a(x) := ^f-(U) is of class 
C 1 except at a point x = x, and vanishes only at x = (by translation invariance). 

Finally, we state the regularity properties for the convex flux, based on the analysis 
in [55], Section 3. The proof is, once again, an adaptation to the general G case of 
the proof of Proposition 3.3 in [22] . which we omit. 

d 2 f 

Corollary C.7. Under convexity of the velocity flux, -3-* > 0, if the shock 
amplitude |it+ — w_ is sufficiently small then the profile is of class C 2 and U'(x) < 
for all i£l. 
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